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Abstract 

We find the exact spectrum and degeneracies for the Gross-Neveu model in two dimensions. 
This model describes N interacting Majorana fermions; it is asymptotically free, and has 
dynamical mass generation and spontaneous chiral symmetry breaking. We show here that 
the spectrum contains 2 N / 2 kinks for any N. The unusual y/2 in the number of kinks for odd 
N comes from restrictions on the allowed multi-kink states. These kinks are the BPS states 
for a generalized supersymmetry where the conserved current is of dimension N/2; the N = 3 
case is the N~ = 1 supersymmetric sine-Gordon model, for which the spectrum consists of 2^/2 
kinks. We find the exact S matrix for these kinks, and the exact free energy for the model. 
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1 Introduction 



The Gross-Neveu model describes interacting fermions in two dimensions. It has no gauge fields 
or gauge symmetry, yet it exhibits much of the behavior of gauge theories in four dimensions. 
The coupling constant is naively dimensionless, but radiative corrections result in a mass scale 
to the theory. Thus the theory is asymptotically free and strongly interacting in the infrared. A 
(discrete) chiral symmetry is spontaneously broken, which gives the fermions mass. In condensed- 
matter language, the interaction is marginally relevant, and there is no non-trivial fixed point. 
The Gross-Neveu model consists of N Majorana fermions if) 1 , with i = 1 . . . N. The action 

is 

S = j ' d 2 z[^d^ + ftdft (1) 

where repeated indices are summed over. At the critical point g = 0, the fermions are free, with 
tjj 1 a function only of z and ip l a function only of z. At this critical point in two dimensions, 
the fermions have left and right dimensions (1/2,0) for ip and (0, 1/2) for rp, so the coupling g is 
naively dimensionless. However, the beta function for this interaction is non- vanishing [EJ|. For 
g > 0, the trivial free-fermion fixed point is unstable, and a mass scale is generated. We denote 
this mass scale M. In this paper we study N > 3; for N = 2 the model reduces to the well-known 
massless Thirring model (the Luttinger model in the condensed-matter literature), and for N = 1 
it is free. The action (|l|) is invariant under the symmetry group O(N). In a separate paper, we 
will discuss generalizations to other symmetries. 

Despite a huge number of papers discussing various aspects of this model, there remained an 
important unanswered question: what is the spectrum of particles? Just because the fields are 
known does not mean the spectrum is known: there may be kink states or bound states, and in 
some cases there are no one-particle states corresponding to the fields themselves. In particular, 
it was long known that for even N that the spectrum contained 2 N / 2 kink states in the two spinor 
representations of the SO(N) symmetry algebra The arguments which lead to the existence 
of kinks for even N apply equally well to odd N, but it was not clear how many such kinks there 
were when N is odd. It was shown in |J that there is no consistent S matrix for a single spinor 
representation, and by utilizing a generalized supersymmetry, Witten showed that the kinks are 
in a reducible representation of SO(2P + 1). He went on to conjecture that they are in two copies 
of the 2 p -dimensional spinor representation. 

In this paper we will answer this question and at long last complete the computation of the 
mass spectrum. We show that there are y/2 copies of the spinor representation when N is odd: 
the number of kinks for any N is 2 N I 2 . A non-integer number of kinks means that there are 
restrictions on the allowed multi-kink states; we give a precise definition below. These kinks are 
in BPS representations of the generalized supersymmetry. For N = 3, the model is equivalent to 
the M = 1 supersymmetric sine-Gordon model, so our result shows that these unusual BPS kinks 
appear even for ordinary supersymmetry. To confirm our claims we compute the exact S matrix 
and free energy of this model, and find agreement with the known results in the ultraviolet limit. 

2 The symmetries 

To determine the spectrum and S matrices of the Gross-Neveu model, we need to understand the 
symmetry structure of the model in depth. There are actually four different symmetries which 
we utilize in this paper. These are all discussed in Q; this section is basically a review of these 
results. 
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O(N) symmetry 

The action has a global O(N) symmetry ip l — > ipiipi and ^* — ► XJ^ T ^ . The matrix f7 must be an 
element of 0(N) because the Majorana fermions are real. The existence of this global symmetry 
means that the particles of the model must transform in representations of this symmetry. 

Spontaneously-broken chiral symmetry 

The action (|l]) has a Z2 chiral symmetry ip — ► —ip, Tp — > Tp. However, for g > 0, this symmetry 
is spontaneously broken, because the fermion bilinear a = ftip 1 gets an expectation value. This 
expectation value results in a mass for the fermions. Note that this expectation value does not 
break the O(N) symmetry; a continuous symmetry cannot be spontaneously broken in two space- 
time dimensions. Equivalently, the discrete parity symmetry ip(z,z) — * ip(z,z) is spontaneously 
broken. 

Local conserved charges 

When a model possesses an infinite number of local conserved currents transforming non-trivially 
under the Lorentz group, it is said to be integrable. The integrability results in powerful con- 
straints on the S matrix, which will be discussed below. It was shown in Q that the action 
(|l]) has at least one of these conserved currents. When g = 0, the energy-momentum tensor 
T = ip l dtp l obeys dT n = for any integer n. This no longer holds for g 7^ 0. However, dT 2 must 
have dimension 5 and Lorentz spin 3. All operators of this dimension and spin must be a total 
derivative: 

dT 2 = BA + dB 

Thus there is a conserved current with components (T 2 — A, B) in the Gross-Neveu model. This 
already requires that the scattering be factorizable, and it seems very likely that there is a 
conserved current of dimension 2n (and hence a charge of dimension 2n — 1) for all integer n. 

Generalized supersymmetry 

A very striking feature of the Gross-Neveu model shown in || is that it possesses an extra 
conserved current of spin N/2. This conserved current is 

i=%^ff---r. (2) 

Another conserved current J of spin —N/2 follows from replacing ip with ip. To prove that this is 
a conserved current even for g 7^ 0, one shows that all possible contributions to 8 J are themselves 
total derivatives. In an equation, 

8J = dJ 

For N = 2, J vanishes and this conserved current is the axial symmetry of the massless Thirring 
model (Luttinger model). For iV = 4, the Gross-Neveu model decouples into two sine-Gordon 
models, and the spin-2 current J is the difference of the energy-momentum tensors of the two 
theories. 

For odd iV = 2P + 1, the results are much more surprising. For N = 3, the current is of 
dimension 3/2, and generates supersymmetry transformations. This seems odd in a theory of 
fermions, but by bosonizing two of the three fermions, one indeed obtains the N = 1 supersym- 
metric sine-Gordon model For general odd N, this results in a generalized supersymmetry. 
The full current algebra of these currents seems quite tricky. It seems very likely that at g = 
it is the WBp algebra studied in [Q], which indeed involves a current of spin P + 1/2 along with 
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those of even integer spin. This was explicitly worked out for the P = 2 case (when g = 0) in |J. 
This current algebra involves the spin-4 local current as well. These currents remain conserved 
when g ^ 0, just like the super symmetry currents do when P = 1. 

Luckily, to determine the S matrix we do not need the full current algebra: we need only to 
understand how the conserved charges act on the particles. The conserved charge Q of dimension 
N/2 — 1 is defined by Q = J dzJ + f dz J, and likewise for Q. The operator Q 2 must be of 
dimension and Lorentz spin N — 2. For example, for N = 3, using the explicit form of the charge 
gives Q 2 oc Pl and Q = Pr, with Pl and Pr the left and right momenta ||. For general N, 
this means that when Q 2 acts on a particle with energy E and momentum p, 

Q 2 oc(E+p) N - 2 , 

since this is the only combination of energy and momentum with the correct Lorentz properties. 
It is convenient to define the rapidity 9 so that a particle of mass m has energy E = mcosh# 
and momentum p = msinh#. The symmetry algebra acting on a particle for odd N can then be 
written 

{Q,Q} = 2m N - 2 e^ 9 {Q,Q} oc 2m N - 2 e -^ e 

{Q,Q} = 2Z (3) 

The central term Z acts on the states as Z\a) = m N ~ 2 z a \a). The (real) numbers z a vanish at 
g = 0, but do not otherwise. As we will discuss, non-zero z a occur because of the non-trivial 
boundary conditions on the kink states 0. 

3 The kink spectrum of the Gross-Neveu model 

One interesting feature of the Gross-Neveu model is that the spectrum is far more intricate than 
a glance at the action (||) would suggest. For any value of N, the spectrum includes kink states, 
as pointed out in [||, 0]. In fact, for N = 3 and N = 4, the spectrum does not contain anything 
but the kink states. The spectrum includes kink states because of the spontaneously broken Z2 
chiral symmetry. The kink interpolates between the two degenerate vacua; in two-dimensional 
spacetime there are domain walls between regions of the two vacua. The main result of this paper 
is to answer the question: how many kink states are there? 

The result requires first finding how the kink states transform under the symmetries discussed 
in the last section. We will find the minimal number of particles required to transform under 
these symmetries non-trivially. In section 4, we will work out the S matrix of these particles, and 
then in section 5 we will show that these particles and S matrices give the correct free energy. 

Again, we follow the analysis of ]2[ and study the kinks semiclassically. The fermion bilinear 
a = ip 1 ^ 1 gets an expectation value do; the kink is a configuration of a interpolating between 
+0"o and -co. In the semi-classical approach, one treats a as a classical background field, and 
quantizes the fermions in this background. The fermions interact with the background by the 
interaction 

gaf-f 

Thus in the semiclassical limit, this problem is equivalent to quantizing iV Majorana fermions 
in a kink background. A famous result of Jackiw and Rebbi shows that the Dirac equation 
with this background possesses a single normalizable zero-energy mode /o,/o |J- The zero- 
energy solution is real, while the finite-energy solutions are complex, so the zero-energy state 
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is created/annihilated by a real operator b l , while the finite-energy states are annihilated and 
created by a l n and ajj respectively. Thus the semiclassical expansion for tp % and ^ is 

oo 

^ = + Y J [U< + fn< 

n=0 

oo 

¥ = f~0b l + J2[fn< + M 
n=0 

The canonical commutation relations for require that 

{b\b>} = 25 ij (5) 

This is called a Clifford algebra; the combinations M lJ = generate the SO(N) algebra. In 

this paper, when we use O(N) we generally mean the symmetry group, while SO(N) refers to 
the resulting symmetry algebra. 

The presence of these zero modes means that there must be more than one kink state, because 
one can always act with any of the b l and get a different state. In other words, denoting the kink 
states by \a), some matrix elements 

</3|6>) (6) 

are non- vanishing. The kink states therefore form a representation of the Clifford algebra. The 
simplest possibility is that the b l act like the gamma matrices of SO(N). A gamma matrix (^ 2 )rs 
obeys the Clifford algebra (||), and is an SO(N) invariant where i = 1 . . . N transforms in the 
vector representation of SO(N), while r, s = 1 . . .2 mt<yN / 2 ^ transform in the spinor representa- 
tion^). The kinks must therefore be in the spinor representation(s) of SO(N). This is effectively 
charge fractionalization, because the vector representation appears in the tensor product of two 
spinor representations. Thus the original fermions ip l in the vector representation "split" into 
kinks in the spinor representation(s). 

The properties of the spinor representation(s) of SO(N) are different for even and odd N. For 
even N, there is an operator ^ N+1 = ib x b 2 . . . b which commutes with all the SO(N) generators 
M ij but anticommutes with the Clifford algebra generators. Therefore the 2 Ar / 2 -dimensional 
representation of the Clifford algebra is not irreducible under the SO(N) when N is even. It 
decomposes into irreducible representations with j N+1 = ±1. These representations are each 
2 N / 2 ~ 1 -dimensional, and are known as the spinor representations of the algebra SO(N). The 
entire kink spectrum for even N consists of 2 N / 2 kinks in the two spinor representations. It 
was shown in 0] that including these kinks along with the other particles gives the correct free 
energy. This effectively proves that these are all the kink states. 

For odd N = 2P + 1, the spinor representation of dimension 2 P is irreducible. There is only 
one such spinor representation of SO(2P + 1). Thus it might seem that the simplest possibility 
is that there are precisely 2 P kinks in the Gross-Neveu model for odd N. However, Witten shows 
that there must be more ||]. The key is to look at the discrete symmetry ^* — > — ij) 1 — > — 
One can think of this symmetry as being (— 1) , with the caveat that the fermion number F is 
not defined for iV odd: only (—1) F is. This symmetry is part of the O(N) group, but not part 
of the connected SO(N) subgroup when N is odd. Because of the form (Q), b % — ► —b l under this 
symmetry. Thus there is a symmetry operator A = (—1)^ which obeys A 2 = 1 and A6*A = — 6\ 
In other words, A anticommutes with all the Clifford algebra generators, and commutes with the 
SO{2P + 1) generators. Thus the kink states in an irreducible representation of SO{2P + 1) must 
all have the same eigenvalue of A. The fact that the matrix element (|6|) does not vanish means 



(4) 
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that some kink states have eigenvalue +1, while some have eigenvalue —1. Thus there must be 
more kink states than those in the single irreducible representation of SO{2P + 1). Note that 
for even N, we can identify A with j N+1 , and the multiple representations with the two different 
spinor representations. For odd N, there is only one kind of spinor representation, so it must 
appear at least twice. 

So what happens for odd ./V? The simplest possibility is that there are kinks in two spinor 
representations, each of dimension 2 P |||, [jljj . We will see that this is essentially correct, but 
there is a major subtlety. The key to the answer is in the generalized super symmetry. 

It is very useful to first examine the case N = 3, where the answer is known ||, 0, 



Here the Gross-Neveu model is equivalent to the jV = 1 supersymmetric sine-Gordon model [Q] 
at a fixed value of the sine-Gordon coupling (3$g where the model has an extra 50(3) symmetry. 
One might expect that the kink spectrum would generalize the kink and and antikink of the 
ordinary sine-Gordon model to account for the super symmetry. For the M = 1 supersymmetric 
sine-Gordon model, it seems plausible that the kink and antikink each are a boson/fermion 



doublet, making four particles in all [16]. However, this is not the correct spectrum. The reason 
is the central charge Z appearing in the supersymmetry algebra ((||) with N = 3). For the 
supersymmetric sine-Gordon model, Z does not vanish (7j: it is the integral of a total derivative. 
Z acting on the kink states is non-zero because here the boundary conditions on the field at 
positive spatial infinity and negative spatial infinity are different. 

The consequences of a non-vanishing Z are familiar in supersymmetric field theories. Rep- 



resentations where the mass m = \Z\ are called BPS states [17]. BPS representations usually 
have a smaller number of states than those with Z = []7|] and so are often referred to as "re- 
duced multiplets". When the BPS condition m = \Z\ holds, one can find a combination of the 
generators which annihilates all BPS states. For the generalized supersymmetry, the condition is 

m N - 2 \ 

O- ! Ve^Qj|a)=0 (7) 
for all BPS states \a). 

We need now to introduce a precise definition of the "number of kink states" K. We say 
there are K kink states if the number of n-particle kink states at large n depends on n as K n . 
Put another way, the entropy per kink of a gas of non-interacting kinks is hiK. K is usually an 
integer in field theory: e.g. if there are two one-particle states, there are four two-particle states, 
eight three-particle states, and so on. The kink structure need not be so simple, though. For 
example, consider kinks in a potential V((p) = (cj) 2 — l) 2 . The potential has minima at = ±1, 
so there are two one-kink states, having <^>(oo) = ±1 and 4>(— oo) = =pl. However, there are only 
two two-kink states. In fact, there are only two kink states for any n. Thus the number of kink 
states in this example is K = 1. The example of most importance for this paper is kinks in a 
potential like V{4>) = 4> 2 (cp 2 - l) 2 . Here the minima are </> = 0, ±1. If kinks are allowed to only 
interpolate between adjacent minima (i.e. only between and ±1), then the number of ra-kink 
states doubles every time n is increased by 2. Therefore, K = \/2 in this triple-well potential. 

For the BPS representations of the generalized supersymmetry (^), K = This was 

already known for M = 1 supersymmetry in two dimensions |l2|, 14, |ll]], the N=3 case here. The 



BPS representations was initially given for the thermal ((^1,3) perturbation of the tricritical Ising 
model, which has N = 1 supersymmetry as well as a Landau-Ginzburg description in terms of a 
field 4> with three degenerate minima [fl^] . To reconcile these two facts, Zamolodchikov showed 
how supersymmetry acts non-locally on the V2 kinks in this triple-well potential. Somewhat 
more surprisingly, there are ^/2 kinks in the J\f = 1 sine-Gordon model as well. 
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To give explicitly the action of the generalized supersymmetry for odd N, we denote a kink 
with (j>(x = — oo) = r and 4>(x = +00) = s as \rs). In the triple-well potential, r,s = 0, ±, 
and r — s = ±1. The central charge Z acting on a kink of mass m in a BPS representation is 
m N ~ 2 (r 2 — s 2 ). The BPS representation of the generalized supersymmetry for odd iV is non-local: 
it changes the kink configuration all the way to spatial — 00. On a single-particle state of rapidity 
0, 

Q\r0) =ir{me e ) iN - 2)/2 \-rO) Q\0r) = r(me e ) iN - 2)/2 \0 r) 

Q\r0) = ir{me- e ) {N - 2 ^ 2 \ - rO) Q\0r) = -r(me- e )^ N - 2 ^ 2 \0r) 



One can indeed verify that Q and Q acting on this representation satisfy the algebra (|3|). More- 
over, the discrete symmetry A acts non-diagonally as 

A\r s) = I — r — s). 

A indeed anticommutes with the generalized supersymmetry generators, With the interpretation 
of A as (— 1) F , we see that the kinks are neither bosons nor fermions. Since the action of Q 
and Q is non-local, we must be careful to define their action on multiparticle states as well. For 
example on a three-kink state, Q acts as 

Q{\a) ® ® | 7 » = (Q\a)) ® \(3) ® | 7 > + (A|a» ® (Q\0)) ® | 7 ) + (A|a)) (A|/3>) ® (Q| 7 » 

and likewise for Q. Thus when Q acts on a kink, it flips all kinks to the left of it. For N = 3, 



this reduces to the action of supersymmetry discussed in [12, 14]. 

Of course, the kinks in the Gross-Neveu model must still transform as a spinor of SO(N). 
Thus we see that for odd N, the simplest non-trivial kink spectrum consistent with all the 
symmetries of the theory is for the kinks to be spinors of SO(2P + 1) and interpolate between 
the three degenerate wells of this potential. Each kink is labelled by a spinor index 1 . . . 2 P and 
by a pair of vacua. The number of kinks is \f2~ x 2 P = 2 N I 2 . We emphasize that the generalized 
supersymmetry and the discrete symmetry A require that there be more than the spinor: 2 N I 2 is 
the minimum number of kinks satisfying these constraints. In subsequent sections, we will find 
the S matrix for these kinks, and show that it gives the correct free energy, effectively confirming 
the presence of these BPS kinks. 



4 The S matrices 

In this section we work out the S matrix for the kinks in the Gross-Neveu model for odd N = 
2P + 1. 

This S matrix must obey a variety of constraints. These are easiest to write in terms of rapidity 
variables, defined so that the rapidity 9\ of a particle is related to its energy and momentum by 
E = m cosh 6\ , p = m cosh 6\ . Lorentz invariance requires that the two-particle S matrix depend 
only on the difference 9 = 6\ — 62- Any S matrix arising from a unitary field theory must be 
unitary and crossing-symmetric. Moreover, the integrability of the Gross-Neveu model means that 
the multi-particle S matrix must factorize into the product of two-particle ones. The resulting 
constraint is called the Yang-Baxter equation. Finally, the S matrix must obey the bootstrap 
equations. These mean that the S matrix elements of a bound state can be expressed in terms of 
the S matrix of the constituents. This constraint is explained in detail in 0, for example. Bound 
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Figure 1: The SO(N) Dynkin diagrams for even TV" and odd N 

states show up as poles in the 5 matrix in the "physical strip" < Im{6) < tt. If there is such a 
pole at 9 = 8j in the kink-kink 5 matrix, then the kinks of mass m have a bound state at mass 

m 3 - = 2mcosh(0,-/2) (8) 

When N = 3 and N = 4, the kinks are the only states in the spectrum ||]. For N = 3, 
this was confirmed by the computation of the exact free energy of the j\f = 1 supersymmetric 
sine-Gordon model |]l^]. For the latter, it follows because the N = 4 Gross-Neveu model can 
be mapped to two decoupled ordinary sine-Gordon models, each at the 5C/(2)-invariant coupling 
(Psg ~ * ^ m conventional normalization) ||. This special decoupling happens because the 
algebra 50(4) is equivalent to SU(2) x SU(2). At the SU{2) invariant coupling of sine-Gordon, 
there are no bound states in the spectrum: the only particles are the kinks in doublets of each 
SU(2). In the 50(4) language, these are the spinor representations. 



For N > 4, there are states other than the kinks. They have masses [18] 



mj = 2msm(-^-^j (9) 

where m is the mass of the kink and j = 1 . . . int(N — 3) /2. The first of these states corresponds 
to the particle created by ^) % \ it is in the ^-dimensional vector representation of SO(N). The 
other states are bound states of the fermions. One useful fact to note is that each type of particle 
corresponds to a node on the Dynkin diagram for SO(N), as displayed in figure 1. For even N, 
the kinks correspond to the two nodes on the right, while for odd N, the kink is the node on the 
right. The vector representation is the node on the left. 

The exact 5 matrix of the vector particles for any TV" was worked out in [p^| . Using the 
bootstrap gives the 5 matrix for fermion bound states, but not the kinks. The explicit expression 
for the 5 matrix of particles with mass mi (which are in the antisymmetric and singlet represen- 
tations of SO(N)) can be found in ||l9|| . The 5 matrix for the kinks for even iV is worked out in 



[ 20| , |21| (see also y]). Thus to complete this picture we need to work out the scattering of the 
kinks for odd N, including the fact that they have the additional BPS structure discussed in the 
last section. 

The BPS kinks in the Gross-Neveu model are in the 2 p -dimensional spinor representation of 
S0(2P + 1), and they also form a multiplet transforming under the generalized supersymmetry. 
Since the two symmetries commute with each other, the simplest two-particle 5 matrix invariant 
under these symmetries is of the tensor-product form 

S(6) = S spmor (6) ® S B ps(0). (10) 

The matrices S sp i nor (9) and Sbps(6) are respectively the 5 matrices for particles in the spinor 
representation of SO(2P + 1), and for kinks in the triple-well potential. 



S 



The spinor part of the S matrix was found in Q. The fact that it is invariant under SO(2P+l) 
means that it can be written in terms of projection operators. A projection operator V a maps 
the tensor product of two spinor representations onto an irreducible representation labelled by a. 
Here a = . . . P, where a = labels the identity representation, a = 1 the vector representation, 
a = 2 the antisymmetric tensor, and so on up to a = P — 1. The representation with a = P 
is the representation with highest weight 2/ip, where /Up is the highest weight of the spinor 
representation. We do not give the explicit expressions for the projectors because we will not need 
them; they can be written explicitly in terms of the gamma matrices for the algebra SO(2P + 1), 
which in turn can be written as tensor products of Pauli matrices. The spinor S matrix is 

p 

s spi nor(e) = ^2fa(G)r a (e) (11) 

a=0 

The functions f a (Q) are not constrained by the SO(2P + 1) symmetry, but all can be related to 
fp(9) by using the Yang-Baxter equation. The result is || 

(9 + ivrA 
fp - 1 ~ I^A fp 

g + ivrA(2a + l) 

fp - a - 1 ~ 0-«rA(2a + l) fp - a+1 {U) 

where the Yang-Baxter equation does not determine A. 

The function fp(9) must satisfy the unitarity and crossing relations. These do not determine 
the function uniquely, because one can obtain a new solution of these equations by multiplying 
any given solution by a function F(9) obeying F{0) = 1 and F{6) = F(in — 6). This is known 
as the CDD ambiguity. The minimal solution of the unitarity and crossing relations is a solution 
without any poles in the physical strip. The minimal spinor S matrix of SO(2P + 1) has 

A 1 1 



N - 2 2P - 1 



and |g 



-i 



M ' S r(i-*A + A)r«»+i)A-£) <13) 

Notice that this has zeroes at = i7rA(2b + 1) for all b = ... P — 1. These cancel poles arising 
from fll2|), and ensure that there are no poles in the physical strip any of the f a . The following 
integral representation for fp will be useful in the future: 



f P {9) = exp 



u 2cosh[(2P - l)cj/2]sinh(a;) 



(14) 



The minimal S matrix for BPS kinks in a triple-well potential can also be determined by 
imposing the same criteria. In addition, we can utilize the generalized supersymmetry. Because 
the generalized supersymmetry operators Q and Q commute with the Hamiltonian, they must 
also commute with the S matrix. Since we know how they act on multiparticle states, this is 
simple to implement. In fact, some of the work has already been done for us. The S matrix for 
the tricritical Ising model describes the scattering of kinks in a triple well, and is invariant under 



ordinary supersymmetry [12]. The generalized supersymmetry algebra acting on the states (H) 
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Figure 2: Representing kink scattering by four vacua 



is related to the ordinary supersymmetry algebra by making the substitution 6 - 
Thus in order to commute with the generalized supersymmetry, 

S BPS (6) oc Stci((N - 2)6) or S BPS (e) oc S TC i(-(N - 2)6). 



±(N - 2)6. 



Such an S matrix will automatically satisfy the Yang-Baxter equation as well (any solution 
remains a solution under the scaling 6 — > X6). 

It is easiest to label the S matrix elements for kinks by their vacua. A two-kink configuration 
can be labeled by three vacua. As shown in figure 2, a two-particle S-matrix element can be 
labeled by four vacua because only the middle vacuum can change in a collision. Thus this S 
matrix element S^V describes scattering the initial state rst to the final state rs't. For the triple 
well, the labels r, s, s' , t take the values 0, ±1. The elements of Sbps are then 



si 



r-lr+l) 
rr 



s 



(rr) 

r±l rq=l 



(0) 



(0) 



Sr+lr+li®) 



S r _l r _i(6) 



l) p isinh[A6>] 




(15) 



where 



ftr 



The tricritical Ising model S matrix has A 
symmetry fixes the sign, requiring that 



= cos(— r . 
4 

1/4, so in general A must be ±(N — 2)/4. Crossing 



A = (-l) 



p+i 



N 



\P+i 



4A 



We should note that some extra minus signs appear in the crossing relations because the kinks 



are neither bosons nor fermions [21]. This solution of the Yang-Baxter equation was originally 
found in an integrable lattice model, the hard-hexagon model, where the kink vacua correspond 
to the heights in the lattice model. This lattice model is part of the "restricted solid-on-solid" 
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(RSOS) hierarchy [22|. More general RSOS models let the heights/vacua run over more values; 
these arise in the multi- flavor and symplectic Gross-Neveu models to be discussed in [34]. 
The last step is to find the function B(9). Crossing symmetry requires that 



B{9) = B(m-e) 



while unitarity requires that 



B(9)B(- 



sinh [X9 - f ] sinh [f + \0] ' 



One can easily find the minimal solution of these relations, where B(9) has no poles in the 
physical strip < ImO < tt. For the N = 3 case, the J\f = 1 supersymmetric sine-Gordon model, 
there are no bound states: the BPS kinks make up the entire spectrum ]|]. There are no poles 
in the physical strip for N = 3. This result is confirmed by calculating the free energy IIH], a 
computation we will repeat in the next section. 

In general, we need a B(9) with poles, because the Gross-Neveu model has bound states. For 
P > 1, there are particles other than the kink in the spectrum. These have masses given by 
(|9|), with j = 1 . . . P — 1. The fermion is the j = 1 state, and higher j corresponds to bound 
states of the fermion [18]. For odd N = 2P + 1 all fundamental representations except for the 
spinor appear in the tensor product of two spinor representations. Moreover, the poles in ( |i"2| ) 
correspond to bound-state particles in the correct representations of SO(2P + 1). We thus expect 
all states (other than the kink itself) to appear as bound states of two kinks. It follows from (||) 
and (0) that the bound states appear as poles in B(9) at 9j = «7r(l — 2jA) for j = 1 . . . P — 1. 
Crossing symmetry means that there must be poles at 9j = iit2jA as well. 

We find that 



B(9) 



sinh [X9 



+ 



"" 1( f r 1,/2) sinhg 



n 

6=1 



sinh [i( 



1 + in 
' — in + i 



| I* [, 



'iJr[3 + (i-i)E + iss] r2 [k + W-Vm 

iir A\ privfis nnlps nnrTPsnnnr 



- i4bnA)] sinh [±(0 + i4bnA)] 
+ iAbnA)] sinh [±(0 - iAbnA)] 

1 + (21 - + 4^s] £ [g + ~ 4^a] 

2 + (2Z — 1)^- — r + /gA + 4~Sa] 



The 1/ sinh(A0 + in/ 4) gives poles corresponding to bound states with odd j. The finite product 
of sinh functions gives the appropriate poles for even j. The infinite product of gamma functions 
assures that crossing symmetry is obeyed. Another way of writing B(9) is as 



B{9) 



sinh [f - X6] 



int\P/2} s . nh j-i ^ + .^_ _ ginh j-i ^ + _ 2 ^ A ^ 

sinh [±(0 - ivr + i(46 - 2)ttA)] sinh [I( 



Mr[i + (i 



'2A 



47TiA. 



r[(2Z-i)^ + 5fe]r[i + (2i-i; 



_!_ J 2_ 

4A ' 47riA. 



- i(46 



2)ttA)] 



47TiAj 



!=1 



L > 2A ^ 47riA. 



r[(2Z 



47riA. 



r [1 + (2/ - 1; 



1 

1A 



4-n-iA. 



47riA J 



Here the 1/ sinh(A0 — in/ 4) gives poles corresponding to bound states with even j, while the finite 
product of sinh functions gives the appropriate poles for odd j. With a little tedium one can 
show that the two expressions for B(6) are in fact equal. The following integral representation 
will be useful: 



B{9) 



V2 



Vcosh 2X9 



exp 



°° ^ c i(2P- 

.nrs OJ 



2sinh[(P-3/2)u;] + sinh[(P - l/2)u] 
4sinh[w] cosh[(P- l/2)w)] 



(16) 



11 



It is most easily obtained by multiplying the two expressions for B>(9), a standard rewriting of 
the logarithm of the gamma functions in terms of an integral, and then taking the square root. 

We have thus determined the S matrix (|l0| ) for the BPS kinks in the Gross-Neveu model. 

By using the bootstrap, one obtains the S matrix for all the particles, for example reproducing 
the S matrix of fjM for the fermions fy. To apply the bootstrap, we need to utilize a number of 
results from integrable lattice models. In this context, the process of obtaining a new solution of 
the Yang-Baxter equation from an existing one is called fusion |24|. For example, in a model 
with SU(2) symmetry, one can fuse the S matrices for two spin-1/2 particles to get that for a 
spin-1 particle. In the simplest cases, the fusion procedure utilizes the fact at a pole 9 = 9j of 
the S matrix, the residue (9 — 6j)S(9o) becomes a projection operator. This happens only at 
imaginary values of the rapidity, so this does not violate unitarity. For example, from (|l2l), it 
follows that at 9 = 9\ = wr(l — 2A), all the fj(9\) = except for j = 1. Thus S sp i nor becomes 
a projection operator V\ onto the j = 1 bound states. This means that one can think of the 
fermion of rapidity 9 as a bound states of two kinks, one with rapidity 9 + 9±/2 and the other 
with 9 — 9\/2. In general, consider a case where the S matrix of two particles a and (3 has a pole 
at 9 = 9j. The residue of the S matrix at the pole is the sum of projection operators 



lim 



The S matrix of the bound state j from another particle 7 related to the S matrices of the 



constituents by the formula [23 



S (i7) = fe y/WP*\ S(ai) ( e + 0j/2)S^H9 - 0,-/2) \J2 ~^ Va 



(17) 



Note that the matrices in ( |l7|) are not all acting on the same spaces, so this relation is to be 
understood as multiplying the appropriate elements. The upshot is that one can think of the 
particle with mass rrij and rapidity 9 as the bound state of kinks with rapidities 9 + Oj/2 and 
9 - 6j/2. 

The f a given by (|l2| ) determine in which representation of SO(2P + 1) the particles for a given 
j transform. For example, the particles with mass m\ transform in the vector representation, 
so there are 2P + 1 of them. These are the original fermions. The particles with mass mi 
transform in the antisymmetric representation and the singlet, so there are 2P 2 + P + 1 of these. 
In general, the particles with mass mj are in the a-index antisymmetric tensor representations 
with a = j, j — 2, . . .. The reason for the multiple representations of particles at a given mass is 
likely that they form an irreducible representation of the Yangian algebra for SO(2P + l) |l{|. To 
compute the explicit S matrix for all the bound states is a formidable task; for some results, see 



19]. However, by using fusion we can compute the "prefactor" of the S matrix, which is necessary 
for the computation of the free energy in the next section. The prefactor of the S matrix for 
scattering a particle in representation with highest weight \i a from a particle in representation [x 
is defined as the S matrix element multiplying the projector on the representation with highest 
weight n a + Hb- We denote the prefactor as S a b(9), and for example, for kink-kink scattering the 
prefactor is 

Sp P (9) = B(9)f P (9). 

An important check on this S matrix is that the scattering closes. This means that all the 
poles in the physical strip correspond to particles in the spectrum. In fact, zeroes coming from 
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the BPS part of the 5 matrix cancel spurious poles coming from the spinor part: without the 
BPS part, the bootstrap does not close. This was noted in ||, where it was termed a 'violation" 
of the bootstrap for spinor particles §§ . Solving this problem gives another way of seeing that 
the extra BPS structure must be present. To find the zeros coming from fusing Sbps, note that 
8j for odd j obeys sinh(A#j + wr/4) = 0, while for even j it obeys obey sinh(A#j — in/ 4) = 0. 
Thus at X6 = —iir/A, the 5 matrix in (15) projects onto the states 

|-11), \1 - 1), |00) j odd 

while at X9 = iir/A, it projects onto 

|-1, -1), |00>, |11) j even. 





oc sinh 




sinh 


\6- 


i3n 
~_ 




oc sinh 


»-<< 


sinh 


X9 + 


%3tt 
~8~ 



Using the formula (17) yields the bound-state 5 matrix. One finds that fusing Sbps gives an S 
matrix effectively diagonal in the kink labels, and moreover, it is the same for all labels. Thus it 
contributes only an overall factor to the 5 matrix of a bound state j from a kink 7, namely 



j odd 



j even 



This overall factor is the only effect of the restricted-kink structure on the bound states j = 
1 . . . P — 1. The non- integer number of kink states and the funny rules do not affect the bound 
states. However, these extra factors do result in zeroes which cancel spurious poles occurring in 
the fusion of S sp i nor for odd N. 

There is an algebraic explanation for the simplicity of the kink structure in the bound states. 
Conceptually, it is similar to the cases with a Lie-algebra symmetry, but here it is because 
Sbps has a quantum-group symmetry. In our case of interest, the quantum- group algebra is 
U q (SU(2)), with the deformation parameter q = e 2m ^ 3 . This algebra is a deformation of 577(2), 
and many of its representations have properties similar to ordinary SU (2). In particular, the kinks 
interpolating between adjacent minima are the spin-1/2 representations of U q (SU(2)). However, 
when q is a root of unity, various special things happen. In our case of a triple well, what happens 
is that the spin-1 representation is very simple. A spin-1 kink goes from either the vacuum 1 to 
— 1, —1 to 1, or from to 0. Note that whichever vacua one is in, there is only one vacua to go 
to. Hence the number of spin-1 kinks is K = 1. The triple-well structure can be ignored once 
one fuses the kinks: effectively the BPS structure goes away for the bound states. 

To conclude this section, we note that kinks transforming under the generalized supersym- 
metry have appeared before, in the ordinary sine Gordon model at coupling j3g G = 16tr j \2P + 1) 



25j. In fact, the Gross-Neveu model and the sine-Gordon model at these couplings are closely 
related. The masses nij (but not the degeneracies) in the spectrum are identical, with the soliton 
and antisoliton of the sine-Gordon model having the mass m like the kinks here [^] . Even more 
strikingly, the soliton and antisoliton are representations of the generalized supersymmetry (|J), 
as shown in [2j|. The reason for the similarity is that both are related to a certain integrable 
perturbation (by the (l,l;adjoint) operator) of the coset conformal field theories 

SO(N) k x SO(N)i 



SO(N) 



k+l 
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The k = 1 case is the sine-Gordon model at pg G = 16tt/N, while the k — > oo case is the Gross- 
Neveu model [26|. The models for all k should be integrable, and for odd N the S matrix should 
have the form 

S = S q ® S B ps 

where q is related to k. The S matrix S q is the S matrix associated with the spinor representation 
of the quantum group U q (SO(N)), which as far as we know has never been worked out explicitly 
in general. In the k = 1 sine-Gordon case, the generalized supersymmetry is extended to an 
ftf=2 version, with generators Q±,Q±. This means that for k = 1, S g is in fact oc Sbps- For 
k — > oo, q — ► 1 and the quantum-group algebra reduces to the usual Lie algebra SO(N). Thus 

liniq_>i Sq — 5 spinor ■ 

5 The free energy 

In this section we compute the free energy of the Gross-Neveu model at finite temperature using 
a technique called the thermodynamic Bethe ansatz (TBA). This gives a non-trivial check on the 
exact S matrix. In the limit of all masses going to zero, the theorem of p^ ] says that the free 
energy per unit length must behave as 

ttT 2 

limF= —c uv (18) 

m— >0 D 

where cjjy is the central charge of the conformal field theory describing this ultraviolet limit. 
The number cjjv can usually be calculated analytically from the TBA, because in this limit the 
free energy can be expressed as a sum of dilogarithms. The cjjv computed from the TBA must 
of course match the cjjv from the field theory. This provides a very useful check, because the 
spectrum is essentially an infrared property. Finding the central charge exactly in the ultraviolet 
limit requires knowing not only the exact spectrum but the exact S matrix as well: the gas of 
particles at high temperature is strongly interacting. All particles contribute to the free energy, 
so if some piece of the spectrum is missing, or if an S matrix is wrong, the TBA will not give the 
correct cjjv 

The free energy is computed from the TBA by two steps. First one finds how the momenta 
are quantized when periodic boundary conditions are imposed. Precisely, one demands that 
the multi-particle wave function obeys ip(xi,X2, ■ ■ .) remain the same when any coordinate x% is 
shifted to x% + L. In an interacting theory, the quantization involves the S matrix, because as 
one particle is brought around the periodic world, it scatters through the other particles. In the 
continuum, this leads to a constraining relation between the densities of states and the actual 
particle densities. The free energy at temperature T is found in the second step by minimizing 
it subject to the constraint. The detailed procedure for this computation has been discussed in 
many places, so we will not repeat these explanations. Several papers closely related to the current 



computation are |2^, 15, 1C]. The TBA computation here is technically complicated because the 
scattering is not diagonal. This means that, as a particle is going around the periodic interval of 
length L, it can change states as it scatters through the other particles. The way to proceed then 
is well known. First, one has to set up the system of auxiliary Bethe equations to diagonalize 
the transfer matrix, then use its results to determine the allowed rapidities. This computation 
amounts to a standard Bethe ansatz computation. The end result can be written conveniently by 
introducing extra zero-mass "pseudoparticles" or "magnons" to the constraining relations. Then 
one minimizes the free energy to find out the equilibrium distributions at temperature T, and 
thus the thermal properties of the 1+1-dimensional quantum field theory. 
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Two kinds of additional complications occur here. The first difficulty is that, as a striking 
consequence of the bootstrap analysis, particles appear in all j-index antisymmetric tensor rep- 
resentations IIj , but in general there is not a single mass associated with a given representation. 
Instead, particles with mass m a appear in representations j = a, a — 2,..., a mixture we call 
p a . On the other hand, little is known about the diagonalization of transfer matrices acting on 
products of SO(2P + 1) representations. To proceed, we have to assume that the usual Bethe 
equations based on the Dynkin diagram of the underlying Lie algebra |28| apply, in the 50(2^+1) 
case, not to the irreducible representations IT,-, but precisely to the mixtures p a . This is a key 
technical assumption, which is checked a posteriori, because the number of particles of each mass 
can be read off from the TBA equations. 

The second is that the kinks in the Gross Neveu model are not only in the 2 p -dimensional 
spinor representation, they also form a multiplet under the generalized supersymmetry transfor- 
mations. This means that instead of one, one gets two auxiliary problems, one to diagonalize 
the transfer matrices acting on the degrees of freedom transforming in the spinor representation, 
the other to diagonalize the transfer matrices acting on the BPS degrees of freedom. A simi- 
lar but simpler problem has been solved in the TBA calculation of the N = 1 supersymmetric 
sine-Gordon model [15|. 



We first solve the auxiliary problem for the BPS S matrix. This was already done in [29], 
because Sbps is the same as that in the tricritical Ising model, up to a rescaling of 9. The result is 
that the TBA equations in the tricritical Ising model require introducing one pseudoparticle with 
density t(9), and hole density f{9). Then the density of states for the particle in the tricritical 
Ising model is 

1 ^ ^k" ^ 1 

2irP TC i(9) = m cosh 9 + -— — * Ptci + ( r ~ f ) 

Z Z7T Z 

where we have defined convolution as 

/oo 
a(0 - 9')b(9')d9' . 
-oo 

The kernel £(6>) is defined as 

2P — 1 

m - c swp-m (19) 

in general, with P = 1 for the tricritical Ising model. The density of real particles is denoted by 
PTC I- It is most convenient to give most kernels in terms of Fourier transforms, defined as 

SO that £(uj) = 75 T-Tnrr . 

sv / 2cosh(^-J 

The density of states for the pseudoparticles is then related to the density of real particles 
Ptci by 

2vr (r + f) = £ * p TCI 

The pseudoparticles have zero energy, but they contribute entropy to the free energy. The 
equilibrium values ptci an d r are found by minimizing the free energy subject to the above 
constraints. Notice that simple manipulations allow us to reexpress the density of states for the 
particle as 

2itPtci(9) = m cosh 9 H — - — * ptci — £, * f = m cosh 9 + £ * r 

2lT 
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The second form will be the most convenient in what follows. 

Before giving the full answer for the odd- TV" Gross-Neveu model, let us first also review the 
solution for the P = 1 case, the N = 1 supersymmetric sine-Gordon model at its SO(3) symmetric 
point. As discussed in the last section, the S matrix is 

-950(3) = Stci ® S^I^Stt- 

The latter piece is the S matrix for the ordinary sine-Gordon model at an 5C/(2)-invariant point 
(Psg ~ * <^ 7r * n ^ ne usua l normalization). The auxiliary problem for this piece is equivalent to 
diagonalizing the Heisenberg spin chain, a problem solved by Bethe 70 years ago. One must 
introduce an infinite number of pseudoparticles with densities Pj(9), in addition to the pseu- 
doparticle density t{9) required for the BPS S matrix. The density of states Pq for the kinks of 
mass m is then related to the particle density po given by 

oo 

2ttP (9) = m cosh 9 + y * p o (0) - CJ i°° ) * Pj( e ) ~£* f ( d )- ( 20 ) 

3=1 

The kernels crj°°'' follow from the Bethe ansatz analysis, and are 

=e -iH. 

The kernel y comes from two places. There is a contribution from the prefactor of the S matrix, 
and for kink-kink scattering there is an extra piece arising from the Bethe ansatz analysis p9|| . 
The formula for any P is 



y^ = ^B + ±i^f P + \((] 2 ,2.,, 



so that 

y-y PP Ip=i-i 4cosh 2 (w/2) 

The other Bethe equations relate the densities of states for the pseudoparticles to particle and 
pseudoparticle densities. They are 



2n Pj (9) = af°Kp Q {9) - (22) 



1=1 

where the density of string states Pj is 

p j = Pj+Pj- 

The new kernel is given by 

AJ°) = 2coth(a;) e -"w*(j.Okl S inh(mm(j,/)a;). 

These equations can be simplified greatly by inverting the matrix Aji (for details see for example 
11 ), giving 

2nPj (9) = 5 j0 m cosh 9 + £ * (pj^ + pj_ x ) (23) 
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Figure 3: The TBA for the 0(3) Gross-Neveu model (the super symmetric sine-Gordon model). 



Even the pseudoparticle r coming from the BPS piece is included in this equation, by defining 
p_i(#) = t(8) and P-i = r + f. Minimizing the free energy then yields the TBA equations 

€j (0) = 5 j0 m coshO - £ * [ln(l + e~^~ xtT ) + ln(l + e -^+ l/T ) 

with j = —1,0,1,2,..., and e_2 = °°- These TBA equations are conveniently encoded in the 
diagram in figure 3. The circles represent the functions e a ; the filled node represents the fact that 
the equation for eo has a mass term. The free energy per unit length F is 

/OO JQ . . 

— cosh<91n(l + e- eoW/T ) (24) 
-oo 2?r v 1 

One can check that in the limit m — » 0, this free energy yields the correct central charge 3/2 |i~5| ]. 
This confirms the presence of the BPS kinks in the spectrum. 

The TBA for the full 0(2P + 1) Gross-Neveu model is conceptually similar to the 0(3) case, 
but is much more involved technically. To complicate matters further, there are particles with 
masses m a , a = 1 ... P. Luckily, many of the technical complications have already been solved. 
The diagonalization of the BPS part is the same as in the tricritical Ising model [29]. This 



requires introducing the pseudoparticle density t(9), as described above. The diagonalization of 
the auxiliary problem for the SO(2P + 1) spinor part is done using the standard string hypothesis 
based on the Dynkin diagram |3(|. For the TBA for the SU(N) and 0(2P) Gross-Neveu models 
discussed in ]1Q[| , the TBA equations are related to the SU(N) and SO(2P) Dynkin diagrams. 
However, SO(2P + 1) is not simply laced, so this case is somewhat more complicated, obvious. 
Nevertheless, the appropriate computation has been done already in |3(J, utilizing the Yangian 
structure of the S matrix. The densities of particles with mass m a are defined as p a o, while the 
auxiliary problem requires introducing pseudoparticle densities p a j with j = 1 . . . oo, as well as r. 



In our notation, we try to remain consistent with [10], where the TBA for the 0(2P) Gross-Neveu 
model is discussed. 

The quantization conditions for the densities of states P a j = p a j + p a j are then 

P oo 

2vrP aiO (0) = m aC o S hd + Y,y^ +1/2) *Pb,m-Y. a[ r ] *~P^ 9 ) 

6=1 j=l 

P oo 

27rPp i0 = mcoBhe + ^yS +1/2) */%,o(e)-^ff^U/5pj(e)-^*f(fl) (25) 

6=1 j=l 

The last terms in equation ( p5| ) are expressions of the logarithmic derivatives of the transfer matrix 
T eigenvalues. The kernels y a b are simply related with the S matrices described in the previous 
sections. Except for the case a = b = P given above in (21), one has 3^f +1 ^ = \ w ^ n ^ab, where 



S a b is the prefactor for scattering particles belonging to the representations with highest weights 
p a and pb- As explained above, these S matrix elements do not probe the kink structure, so the 
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same computation will give these kernels for O(N) for even and odd N. Thus these kernels for 
odd N can be read off from the N even results in |1Q|| , giving 

1^1 cosh((P — 1/2 — max(a, b))uj) sinh(mm(a, b)u) 



y, 



(P+l/2) 
ab 



UJ 



Sab 



y a,P 

f;(^+l/2) 
Spp 



cosh((P- l/2)u>) sinh(u;) 
sinh(au;) 



1 



2cosh((P - l/2)w) sinh(a;) 
\u\ sinh(Pu;) 



4cosh((P - l/2)w) sinh(a;) cosh(w/2) ^ 

for a, b = 1...P — 1. The latter kernel came from (^). The pseudoparticle densities are given 
by the auxiliary Bethe system 

oo P—l oo 
2vr p a j(9) = crj°° ) * Pa,0 *K ab *p b ,i{0) -^A}ff 2 *KaP*pP,l 



2irp PJ (9) 



2tt (t + f) 



(oo) 

v) /2 J *pp,o 

£ * PPfi 



1=1 b=l 
oo P-l 



i=l 6=1 



Z=l 

oo 

_ v 4 (oo) 

j/2,Z/2 

i=l 



* -Kpp * pp 7 i 



(27) 



where the kernel K in Fourier space is 

K a , a ±i = 
Kpp-i = 



1 



2 cosh 
coshu;/2 



Kp P 



1 

cotha;/2 
cothu 



(28) 



cosh uj 

The equations (|27| ) are the continuum limit of the SO(2P + 1) Bethe equations, with sources 
terms associated to all the p a representations. In the preceding equations (p^-28), the index a 
takes values 1 ... P — 1; the equations involving the index P are given explicitly. The key feature 
of this system is the appearance of factors of 2 in the terms involving the latter, which in the 
Dynkin diagram corresponds to the spinor representation. It is directly related to the fact that 
the P th root is the shortest, and has length 1 while all the others have length 2. 
The TBA equations are written in terms of functions e a j, defined as 

Pa,j _ 1 t _ 1 

P a ,j ~ 1 + e e ^il T t + f ~ 1 + e ep --^ T ' 

for a = 1 ... P. The values j runs over depends on the value of a: for a = 1 ... P — 1, j takes 
values 0, 1, . . . oo, while for a = P, j = —1, 0, 1 . . . oo. The extra function ep_i arises from the 
diagonalization of the BPS S matrix. This problem can now be put in a considerably simpler 
and rather universal form by inverting the kernels. The resulting TBA equations are 

-T(f> * [In (1 + e - t ^+ 1 / T ) + In (l + e ~ e ^- l/T ^ 



together with 
ep-i,#) 



+T(p * [In (1 + e e "+^/ T ) + In (l + e ta -^ /T 



1, 



P-2 



(29) 



+TV>*m(l + e ep ^ /T 



In (l + er^-i.i+iAj + In (l + e -*p-w-i/ T 
In (l + e e r-*,i/ T \ + In (l + e ep > 2 J+ l/T ) + In (l + e ep ^- l/T 



(30) 
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Figure 4: The TBA for the 0(2P + 1) Gross-Neveu model (here P = 4) 



and 

e Pj -(0) = -Tf* In (^1 + e" ep ^+ l/T ) + In (^1 + e - £p ^- l/T ) + T£ * In f 1 + e ep - x ' i >^ T '\ (31) 

In the last equation, the coupling to tp-i.j/2 occurs only when j is even. The kernels are defined 
by their Fourier transform, 

I 1 7 cosh(f) £ 1 



2cos1i(lj) cosh(o;) 2cosh(^) 

As usual, the mass terms disappear from the equations, but are encoded in the asymptotic 
boundary conditions 

£a,o - ► "T-a cosh ep f o — ► m cosh 9, 6 — > oo 

The TBA equations can be conveniently encoded in the diagram of figure 4. The is the non- 
simply-laced generalization of the diagrams of fTOfl . One remarkable feature is how the extra 
pseudoparticle coming from the BPS structure (the node on the bottom left) fits in perfectly 
with the pseudoparticles coming from the diagonalization of the spinor part of the S matrix. In 
this sense one could infer the existence of the BPS structure from a careful examination of the 
structure of the TBA equations. 

The entire purpose of this section is to verify that this S matrix does give the correct free 
energy ( [1S| ) in the ultraviolet m — > 0. This is of course a major check: it shows that this 
complicated spectrum and S matrix conspire to give the correct central charge c = -P + l/2in 
the UV, corresponding to 2P + 1 Majorana fermions. The free energy per unit length F is given 
in terms of these dressed energies e a as 

— cosh0 1n(l + e- ea ' o(e)/T J (32) 

-oo 27T V / 



a 



By rewriting F as m — * in terms of dilogarithms, we have verified that one does indeed obtain 
c = P + 1/2 as required. Thus we take it as proven that this is the correct S matrix. 



19 



Another check is that in the IR limit m a — > oo, the correct particle multiplicities are obtained. 
Namely, one should have 



lim F(m,T) = —T'S^ n„m a [ — cosh#i 



e -nia cosh(0)/T ^23^) 



where n a is the number of particles with mass m a . To obtain this by taking the limit of ( |32| ) 
takes a little bit of work. Define 

Y aj = lim e £ <^ (e) j > 1 

Y a0 = lim e ««,oW-m„ooBhe 

It then follows from ( ^J3^ ) that n a = Y a ,o- To find the Y a j requires taking the m — > oo limit 
of the TBA equations (|2^,30|,31). In this limit, one can replace the e a j(6) by ]nY a j. Thus the 
integrals can be done explicitly, giving a set of polynomial equations for the Y a> j. For example, 
for j = they are 

(Y a , ) 2 = Y a _ li0 Y a+li0 (l + Y aA ) a = l,...,P-2 

(yp-i,o) 2 = i r p-2 I o(y P ,o) 2 (i + ^p-i,i)yp,i(i + ip,i) -1 ip,-i(i + >p,-i) _1 

(y P , ) 2 = F P _ lj0 (i + Fp,-i)(i + yp,i) 

where Yq o = 1. The (not written) equations for Y a j for j > 1 do not depend on the 3^o> 
so these can be solved separately. One finds immediately that Yp—\ = 1, but unfortunately, 
we were not able to derive an explicit closed-form solution for the other Y a j. However, it is 
easy to find them by solving the polynomial equations numerically. For example, for P = 2, 
one finds that Y\^ = 14/11, and Yi t \ = 11/5. Plugging this into the equation for Y a ^ gives 
^i,o = n i = 5 and Y^o = ri2 = 4\/2- These indeed are the correct multiplicities for P = 2. For 
P = 3, we find Y 1A = 27/22, Y 2 ,i = 95/147 and F 3j i = 21/11. This yields n x = 7, n 2 = 22 
and n3 = 8\/2- In particular, note that ri2 = 21 + 1, the dimensions of the antisymmetric and 
singlet representations of 0(7). This checks that the transfer matrix diagonalization in |?0| indeed 
includes all the particles at each mass, even though they come from more than one representation 
for o = 2...P-l. 

As another check, we can easily generalize this TBA calculation to the SO{2P+\)j t x SO(2P+ 
l)i/SO(2P + l)fe+i coset models mentioned at the end of the last section. The TBA is almost 
the same as that above, except that the right hand side is truncated, so the j in e a j runs only 
from 0. . . k — 1 for a = 1, . . . ,P — 1, while j = —1, 0, 1 . . . 2k — 1 for the a = P nodes |3(|. One 
can check that for this truncated TBA system 

k(2P + 1)(4P - 1 + Jfe) 



cuv 



2(k + 2P- l)(Jfe + 2P) 



as required. For k = 1, the conformal theory has central charge c = 1, and the coset perturbation 

3 2 

coincides with the sine-Gordon model for -# 

07 

in figure 5, and was first studied in p5l Ej. 



02 

coincides with the sine-Gordon model for = 2 p +1 . This TBA is represented by the diagram 



6 Conclusion 

We have completed the solution of the Gross-Neveu model for any number N of fermions. For 
odd N, generalized supersymmetry results in the existence of BPS kinks. We found the exact S 
matrix for these kinks, and used this to compute the exact free energy. 



20 



Figure 5: The TBA for the sine-Gordon model at ^ = ^rpj (here P = 4) 

One striking feature is that there are a non-integer number K = 2 N ^ 2 of these kinks, in 
the sense that the number of n-kink states goes as K n . In the simplest case, the generalized 
supersymmetry reduces to M = 1 supersymmetry in two dimensions. We note that a number 
of results concerning BPS kinks in 1+1 dimensional supersymmetric field theories have been 
derived recently (see |3!| and references therein), but to our knowledge, none of these papers 
have discussed this multiplet with the y/2 particles. However, unusual particle statistics related 
to Clifford algebras have been discussed in [|33|], 

We have generalized these results to a large number of models with four-fermion interactions 



[34]. These models include Gross-Neveu-like models with Sp(2N) symmetry, and multi-flavor 
generalizations of the Gross-Neveu model. As in the 0(N) case, kinks with non-integer K appear 
whenever the symmetry algebra is not simply laced. 
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